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ABSTRACT 



The amplitude distributions of the sky and sea background noise in the 
3—5.6 fim and 8—14 fim spectral regions are analyzed. These background scenes were 
obtained with an AGEMA 780 Thermovision sensor and recorded with the AGEMA 
TIC— 8000 data acquisition system. The amplitude distributions of 3—5.6 fjxn clear 
sky, 3—5.6 nm cloudy sky, and 8-14 /xm sea bzickground noise sources obey Gaussian 
distributions. A Gaussian amplitude distribution does not provide a good fit to the 
observed 3—5.6 fim sea background data and the 8—14 fim clear sky beickground. 
Scene contamination may have caused the 3—5.6 fjm sea beickgroimd to deviate from 
a Gaussiein distribution. In addition, the eunplitudes of the 8—14 /xm clear sky 
beickground fall into the non-linear portion of the dynamic range of the scanner. As 
a result, no conclusion is made regarding the parent distributions for 3—5.6 /xm sea 
background or the 8—14 /xm clear sky background. 
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I. INTRODUCTION 



A. GENERAL 

During recent yeeixs, a large number of sophisticated Infrared detection 
systems have been investigated euid developed. In an Infrared detection system 
it is necessary to suppress background noise to detect the target. 
The suppression of infrared backgrounds for the purpose of target detection 
requires a knowledge of the statistical metrics of the amplitude distribution. 
Specifically the mean, variance, and parent distribution must be known for 
each type of background. 

In 1974, Itakura et al. [Ref. 1] presented the results of a statistical analysis of 
IR backgrounds obtained with a radiometer over the wavelength range from 2 
to 14 /im. Statistical properties of sky—, forest— and cities— backgrounds were 
described in this study. The background noise was described as a set of random 
two-dimensional pulses. The results showed the amplitude distribution of 
background spatial radiance, for thermal radiation emitted between 8—14 fim, 
obeyed Gaussian statistics whereas in the spectral region dominated by scattered 
sunlight (2—3 ;im) the amplitude distribution followed Poisson statistics. 
The pulse width was described by a Poisson distribution. 

Manolopoulos [Ref. 2 ] measured the radiance of several objects with an 
AGEMA THV 780 infrared scanner. The purpose of this study was to provide 
background radiance values for the development of the AN/SAR— 8 Infrared Search 
and Target Designation (IRSTD) system. In order to validate the Schwartz— Hon 
(SH) model of the sea surface emissivity and the sky radiance incident (or reflected) 
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on the aea, thermal images of the sky and the sea background radiance were 
analyzed by Psihogios [Ref. 3] with the same thermal imaging system. Although 
both of these investigations used the AGEMA THV 780, they did not perform an 
indepth statistical aneJysis of beickground amplitude distributions. 

B. PURPOSE AND ORGANIZATION OF THESIS 

The purpose of this study is to describe the amplitude distribution of 
backgrounds which were not studied by Itakura et al. [Ref. 1]. The backgrounds 
investigaited here include clear sky (3-6.6 fan and 8—14 fan), overcast sky 
(3—5.6 /im), and sea surface (3-6.6 fan and 8-14 /rni). No attempt was made to 
investigate the spatial frequency of background noise. 

The present work consists of five chapters including the introduction. 
The thesis is organized as follows: 

* Chapter I is the Introduction. 

* Chapter II is an overview of the subjects considered during this study such as 
thermal radiation laws and statistical methods. 

* Chapter III describes the equipment and the techniques of measurement. 

* Chapter IV presents a description of the measurements zind the results of 
the data zineilysis. 

* Chapter V discusses the results of chapter IV 2 ind presents recommendations 
for further studies in this area. 

* Appendixes A and B show the listing of FORTRAN programs SBl aind SB2, 
respectively. 

* Appendix C describes the transfer procedure of the thermal image data 
collected using zin if800 microcomputer. 
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n. BACKGROUND TO PROBLEM 



A. IR RADIATION AND THERMAL IMAGING 

All non-target objects which scatter, reflect, or emit infrared radiation are 
considered background noise sources. Examples include the sun, clouds, land, sea, 
sky and birds. Measurement of the background noise can be done using 
a "Thermal Imaging System". Thermal imaging systems are thoroughly described 
by Lloyd [Ref. 4]. 

Four wavelength regions of the IR spectrum can be used in IR imaging. These 
include the near infrared region with wavelengths of about 1 fan and other windows 
at 2—2.5 fan, 3—5.6 fan, and 8—14 fan. Detection in other parts of the spectrum is 
restricted by atmospheric absorption by water and carbon dioxide molecules 
(see Figure 1). Imaging in all these regions can be carried out with ambient or 
active illumination. For wavelengths greater than 3 fan, the object's black body 
thermal radiation is sufficient for thermal imaging. However, the contrast of 
thermal imaging is very low, particularly in the 8—12 fan portion of the spectrum 
where most of the thermal radiation energy is concentrated. Since a thermal 
resolution of 0.1 K is typically sought with background at 300 K, the uniformity 
requirements of image sensor responsivity in this region are very severe. [Ref. 5]. 

The "AG A 780 Thermovision", which was used to collect the data, is a Uquid 
nitrogen cooled thermographic device able to detect IR radiation in the Middle IR, 
(3 — 5.6 fan designated SW by AGEMA), and in the Far IR, (8 — 14 fan 
designated LW) regions. Although these beindwidths do not coincide with 
the atmospheric window shown in Fig. 1, the overlap is considered sufficient. 
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Figure 1. Spectral Transmittance of the Atmosphere (Reproduced from Ref. 1) 

B. THERMAL RADIATION LAWS 

In 1860 Kirchhoff introduced the law that states, in effect, that good absorbers 
are also good radiators. He also proposed the term blackbody to describe a body that 
absorbs all of the incident radiant energy and that, as a consequence of his law, 
must also be the most efficient radiator. In 1900 Planck described the spectred 
distribution of the radiation from a bleickbody as 



where is the spectral emittance, W cm'^/rm'^ A is wavelength, fan, T is the 
absolute temperature, K, Cj is the first radiation constant equal to 
3.7415xl0^±0.0003 W cm'^fon^ and is the second radiation constant equal to 
1.38054*0.00018x10-'^ W sec K'K [Ref. 6]. 




( 2 - 1 ) 
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Planck's law cam also be written in terms of photon cm”^ sec”* /zm”* by 
dividing Eqn. 2 — 1 by hc/X, which is the energy associated with one photon, giving 



Q = 



A4 -1) 



( 2 - 2 ) 



where Q is the spectral radiant photon emittance, photon sec”* cm”^, and c} is a 
constant equal to 1.88365x10^ sec'* cm'^ /im^. [Ref 6]. 

Figure 2, which is reproduced from Ref. 6, shows the spectraJ radiant 
emittance of a blackbody at veirious temperatures, ranging from 500®K to 900°K. 




Wavelength (microns) 



Figure 2. Spectral Radiant Emittance of a Blackbody 
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G. THE EFFECT OF BACKGROUND NOISE 

The effect of the sea and sky background noise can be seen easily in Fig. 3. 
Figure 3(a) and Figure 3(b) are surface plots of the long wavelength image of an 
appro 2 u:hing ship. In Fig. 3(a), the peak created by the ship is hardly recognizable 
on the intersection line of the sea sind the sky background and the signal— to-^ioise 
ratio is obviously small. In Fig. 3(b) the target is closer and easily distinguishable; 
the signal— to— noise ratio is also improved. 

It is seen clearly that the ampUtude of the sky background noise is higher than 
the sea background noise because the temperature of the sky is higher than the 
temperature of the sea. Whenever the temperature of either sea or sky changes, the 
amplitude of the noise fluctuations will also change without changing the 
distribution function of the noise. 

D. THE STATISTICAL MODEL OF THE BACKGROUND NOISE 

The statisticeJ properties of infrared background noise have been analyzed 
previously [Ref. 1]. The statistical model developed was a rzuidom set of 
two-dimensional pulses whose amplitudes obey either Gaussi 2 m distributions in the 
8—14 fim region or Poisson distributions in the 2—3 fjcm region. 

It is assumed that the radiance received by the system obeys Gaussian 
statistics in both the 3—5.6 fim and 8—14 fim region, as follows; 

p^(„) = ^ (2 - 3) 

where Pg is a Gaussiein distribution; n is the radiance of a certain point, h is the 
average rzuliance, cr^ is the variance of n. 
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Figure 3. Surface Plot of Apparent Thermal Values (lU): (a) target is distemt; 
(b) Target is close . 
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E CHI-SQUARE GOODNESS OF FIT TEST 



1. Procedure 

In this investigation it is assumed that the background noise comes from 
a Gaussian distribution. The null hypothesis states that there is no difference 
between the expected and actual sample distributions. The null hypothesis is tested 
using a Chi— squared test. After obtaining each infrared image, the data are grouped 
into amplitude cells. The frequency— of— occurrence values have been compared with 
those expected under the null hypothesis. The test procedure uses statistics given 






(observed — estimated expected) ^ 
estimated expected ‘ 



(2-4) 



This expression has approximately a chi-squaure distribution with k-l—m degrees 
of freedom, whereby k is the number of cells, m is the number of parameters 
estimated [Ref. 7]. The test procedure used in this study is: 

* if ^ j • accept null hypothesis 

* < X^ < X^ , withhold judgment 
•Ol.k-l-m ■ -OOS.k-l-iB 

* , reject null h 5 pothe 8 is 

• 006)k-l-iti 

2. Goodness of Fit for Normal Distribution 

If the null hypothesis states that the underlying distribution is normal, 
the X test must be preceded by a choice of cell size and an estimation of the mean 
value and standard deviation. These are determined using 



fi — n , (2 — 5) 



where nj is the raidiance of point j, n is the average radiance and N is the sample 
size. The number of parameters estimated is two and hence m=2. [Ref. 7]. 
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III. EQUIPMENT 



A. THE AGA THERMOVISION SYSTEM 

In this study an AGA 780 Thermovision Thermal Imaging system was used as 
measurement eqxiipment. The Thermovision 780 system, a fourth generation 
development of "AGA Thermovision" which was first introduced in 1965, is 
designed for both image recording and image analysis. The basic Thermovision 780 
system [Fig. 4] consists of a Black/White monitor chassis together with the dual 
scanner unit [Ref. 8]. In addition to this basic system, the Thermovision 780 System 
used in this study also includes a Sony RGB color monitor and an IBM AT 
microcomputer configured with AGEMA system cards and software. 
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B. BASIC DESCRIPTION OF THE SYSTEM UNITS 
1. The Dual Scanner Unit 

This unit consists of two seried sczuining cameras each containing a single 
detector. One camera covers the 3 to 5.6 fim spectral band (SW) and the other 
covers the 8 to 14 /im spectral band (LW). The cameras convert electromagnetic 
thermed energy radiated from an object into electronic video signals which are 
amplified and treinsmitted via an interconnecting cable to the Black/White monitor 
chassis. [Ref. 7]. 

а. Electro— Optical Scanning Mechanism 

The two bands of the scanner use separate optics, which make it 
possible to cover both the long and short wavelength simultaneously. Two rotating 
refractive prisms with eight sides, one on the vertical axis and one on the horizontal 
axis, axe used to scan the scene in a raster pattern. The synchronized rotation of 
horizonted and vertical prisms provides 100 horizontal scanning lines, with four to 
one interlace producing 400 lines per frame. The lenses of the system axe maide of 
germanium for LW aind silicon for SW with 7° fields of view. The lenses can be 
focused in the same way as normal photographic lenses. The aperture stop can be 
adjusted between f/1,8 2 ind f/20 amd is controlled from the front of the scanner unit. 
Another control on the front of the scanner allows the selection of a specific filter. 
The data has been collected using lenses with a 7° field of view, the f/1.8 aperture 
stop and the unfiltered position. [Ref. 8]. 

б. Infrared Detectors and Cooling 

An Indium— Antimonide (InSb) detector in the 3 to 5.6 pm band 
aind a Mercury— Cadmium— Tellur ide (HgCdTe) detector in the 8 to 14 pm band axe 
used. Both axe single element point detectors aind are cooled to — 196°C using liquid 
nitrogen. [Ref. 8]. 
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c. Technical Data 



Technical data for the dual AGA 780 scanner is reproduced from 
Ref. 8 and Ref. 9 and shown in Table 1. 



Table 1. AGA 780 TECHNICAL DATA 



Performance Characteristics 


Spectral Region fgm) 


5.6 and/or 8-14 


Spectral rate (Sec"*) 


6.25 


Field rate (Sec"*) 


25 


Interlace 


4 to 1 


Replaceable Fore-Optics 


FOV (deg.) AZ 


7.0 


EL 


7.0 


IFOV (mrad.) AZ 


1.1 


El 


1.1 


NEAT at 220C (^C) 


0.12 


MDT 


< 0.1 


Dynamic Range (OC) 


- 20 to 900 


Optical Data 


Effective Aperture Area (cm^) 


: 24.0 


Diameter of Aperture fcm) 


5.5 


Effective Focal Length (cm) 


: 9.9 


f/ number 


: 1.8 


Detector and Cooler Data 


Type of Material 


InSb or HgCdTe 


Number of Elements 


1 


Peak Wavelength 


5 or 10 


Cooler Type 


liquid nitrogen 



2. The Black/White Monitor Chassis 

The video signal from the scanner imit is amplified and processed within 
the monitor chassis and applied to the display screen. The monitor chassis also 
supplies the required sweep signals, isotherms, and power supply. Output from this 
unit IS connected to the IBM AT microcomputer via appropriate interface cable. 
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The Black/White monitor screen, which displays the 280 lines of 400 
lines per frame, is covered by a mask [Fig. 5] with vertical and horizontal scales 
which are illuminated by the display raster. The horizontal numbers at the top of 
the screen show the Thermal Range setting, from 2 to 1000. The vertical scale, 
graduated —0.5 to +0.5, corresponds to the fraction of the range selected. [Ref. 8]. 

Table 2 describes the operation of the controls associated with the 
conversion from isotherm iinits to radiance and Figure 6, reproduced from Ref. 8, 
shows the controls and indicators on the front panel of B/W monitor chassis. 




Figure 5. B/W Monitor Screen Mask 
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Table 2. BLACK/WHITE MONITOR CHASSIS 
SELECTED FRONT PANEL CONTROLS 



Item 


Description 


Function 


1 


POWHl 


This push UN/UFF button controls the 
ac supply to the THV 780 system. 


9 


THERMAL RANGE 


This nine position switch selects the 
thermal span of interest. The switch 
is calibrated between 2 and 100 lU. 


10 


THERMAL LEVEL 


This control sets the thermal level of 
the thermal image. 


11 


THERMAL L. ADJ 


This control provides fine adjusment 
for calibration of the thermal level 




Figure 6. Black/White Monitor Chassis Controls and Indiceitors 
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3. The IBM AT Microcomputer 



An IBM AT computer is configured with AGEMA system cards and 
softwcire in order to build the Thermal Image Computer (TIC— 8000). This 
TIC— 8000 has menu-directed AGEMA software and works with the MS-DOS 
operating system [Ref. 10]. Previously, the BMC if800 microcomputer, DISCON 
scan converter and IR— Link had been used for the same purposes. The TIC— 8000 
m 2 ikes possible live mode image display 2 uid processing as well as storing the data 
into binary files which can be euicessed easily. The TIC— 8000 with RGB monitor 
and Cannon color printer has been used for the initizd analysis and digital recording 
of data. Unfortunately, there exists a problem in the current configuration. RANGE 
and LEVEL do not transfer to the TIC— 8000 automatically. They must be entered 
by hand for each stored image. 

а. Digital Image 

The amplified analog video signal from the B/W monitor chassis is 
converted to an 8 bit digital signal with the full binary range of 256. The 140x140 
pixel element matrix from a single vertical scan is processed in real time and 
displayed on the RGB monitor [Fig. 7]. 

б. Data File Structure 

The image czui be stored into the image sub directory as a bineiry 
file with a record length 20446 bytes where the image data starts at byte 846. The 
140x140 pixel element matrix is stored in such a way that the even numbered lines 
are stored in the first 9800 bytes and the odd numbered lines in the second 9800 
bytes. 
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Figure 7. The TIC— 8000 Screen with The 140x140 Pixel Digital Image 
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C. MEASUREMENT TECHNIQUES 

The AGA THV 780 infraired scajiner measures the scene reidiance eind stores 
the result in Isotherm Units (lU) which is a arbitrziry unit of measurement. 
The detected radismce for a given pixel is assigned a thermal value. The thermal 
veilue is proportion 2 il to the received photon radiation and is related to the 
blackbody temperature through a nonlinear curve defined by individual calibration 
constants. [Ref. 8]. 

The two basic methods of measuring the temperature by using AGA THV 780, 
direct measurement and relative measurement, are thoroughly described in Ref. 8. 
The direct measurement technique has been used to measure the scene radiance in 
this study. 

1. Direct Measurement 

The equipment is used to obtain the scene radiance directly without the 
use of an external reference source of radiation. The thermal value is set by the 
values of the "Thermal Level" and "Thermal Range" knobs on the B/W monitor 
chassis where the thermal range sets the dynamic range of the scanner and the 
thermal level sets the median value of thermal value to which the scanner responds. 
After adjusting the "Thermal Level" and "Thermal Range" controls to get a 
satisfactory picture end the "Isotherm Level 1" control to highlight the point of 
interest, the marker reading of isotherm scale is multiplied by the "Thermal Range" 
setting and added to the "Thermd Level" setting. The result is the measured 
thermal value in Isotherm Units (lU) which represents the scene radieince [Fig. 8] or 
alternatively the scene rawdiant emittance. Assuming the object is a blaickbody and is 
positioned at a distance of one meter from the scanner, the calibration function can 
be used directly to find temperature [Ref 8], 
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2. CalibmtioD Fxinctions 

The calibration cvirves rised to convert measured thermal value to object 
tempereiture are described by the mathematical model: 



I = 



A 

C 



(3-1) 



Here I is the thermal value, T is the absolute temperature (K) aind A,B,C are 
factory calibration constants. The individual calibration data is given in Table 3 
and the calibration curves are illustrated in Figure 9. [Ref. 8]. 
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Tftblc 3. INDIVIDUAL CALIBRATION OF 780 DUAL SCANNER 
(Reproduced from Ref. 8) 



SERIAL NUMBERS 


SCANNER 


SW» 4011 


LWDB 4011 


DETECTOR 


G 1739 


G 1739 


FILTER 


NOF 


NOF 


LE2^S 


7 3009 


7 3105 


CALIBRATION CQNDniONS 


AMBIENT TEMPERATURE 


23 °C 


24 oC 


RELATIVE HUMIDITY 


50 X 


65 X 


OBJECTIVE DISTANCE 


1.0 m 


1.0 m 


CALIBRATION DATE 


82-09-17 


82-09-27 


OPERATOR 


B.A 


B.A 


CALIBRATION 


CURVE CONSTANTS 




A 


811203 


-3581 


B 


3117.74 


1506.49 


G 


1.00 


-0.436 
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Because thermal veilue is a non-standard quantity, the images sire 
converted to values of radiant emittance. There exists a linear relationship between 
thermal value eind radiant emittance [Ref. 5]. This relationship can be written as 

W = m.[I + c] (3-2) 

where W is radiant emittance (W cm"^), I is thermal value (lU), and m and c are 
mdi\ddual calibration factors. 

The reidiant emittance can be calculated directly from the Eq. 3—2 once 
m and c have been determined. In order to determine these constants, a FORTRAN 
program called SBl has been written (listing of this program appears in 
Appendix A). The program first calculates the temperature from the thermal vedue 
using equation 3—1. These temperatures are then used to calculate radiant 
emittance by integrating Planck' Law over the specified spectral region. The results 
of SBl are shown in Figures 10 and 11 for S\V and LW, respectively. The linear 
relationship between radiant emittance and thermal value can easily be observed in 
both Figures. By usmg a simple Imear computation, these relationships were 
defined as 

\V = 4.69x10-5 • [I + 1] for SW (3 - 3) 

and 

W = 3.03xl0-« •[! + 3.624] for LW. (3 - 4) 

Thus, the parameters of Eq. (3 — 2), m and c, are respectively 4.69xl0"5 (W cm"^ 
IU-‘) and 1 (lU) for SW, 3.03x10"' (W cm'2 lU'^) and 3.624 (lU) for LW 

These calibration constants assume an object distance of one meter. 
Because the distances of the background noise sources are poorly defined, an 
apparent radiant emittance is introduced. The apparent radiant emittance is defined 
here as the radiant emittance of a blackbody at a distance of one meter 
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Figure 10. SBl Output for SW (4 lU to 50 lU) 




THERMAL VALUE (U) 



Figure 11. SBl Output for LW (19 lU to 70 lU) 
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IV. MEASUREMENTS AND DATA ANALYSIS 



A. MEASUREMENTS 

The measurements were mewie from the roof of Building 232 (Spanagel Hall) of 
the Naval Postgraduate School and from Hopkins Marine Station. The 
measurements from the roof of Speinagel Hall were made between 12 and 4 o’clock 
P.M.. In addition, images of the LW sea surfzice were taken during Nov, 1988, from 
the Moss Landing Marine Lab [Ref. 3]. These images were transferred from the ifSOO 
computer to the IBM— AT computer using the DISPORT transfer program (the 
description of the image transfer procedure appears in Appendix C.). Table 4 shows 
the list of measured objects with initial measurement parameters. 



Table 4. LIST OF MEASURED OBJECTS 



Date 


Band 


# of Img. 


Object 


Air Tmp. 


Veather 


Az.El. 


ii.ol.dd 


LW 


30 


Sea 




— 


Random 


11.05.88 


LV 


45 


Sea 


19 OC 


— 


Random 


08.31.89 


SW 


19 


Sky 


20 °C 


Clear 


Random 


09.01.89 


SW 


08 


Sky 


19 OC 


Cloudy 


Random 


09.05.89 


SW 


10 


Sl:^ 


20 t>C 


Clear 


SE,450 


09.05.89 


SW 


05 


Sky 


20 t>C 


Clear 


S ,700 


09.05.89 


SW 


05 


Sl^ 


20 OC 


Clear 


SV,750 


09.08.89 


LV 


20 


Sky 


19 C 


Clear 


Random 


09.14.89 


SV 


20 


Sea 


18 C 


— 


Random 


09.15.89 


LV 


20 


Sky 


21 C 


— 


Random 
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B. DATA ANALYSIS 



1. Procedure 

Because of the vignetting ait edges of the field— of— view, during the data 
analysis only 40x40 pixels from the central part of the images were used. This 
reduced the FOV to 2° and the instrument factor to a negligible level. This also 
reduced the sample size from 19600 to 1600. 

Each sample value is converted to apparent radiant emittance (W cm~^) 
from Isotherm Units (lU). A FORTRAN program called SB2 was written for this 
conversion (listing of this program appears in Appendix B). SB2 reads the image 
file and writes the output as apparent radiant emittance (W cm~^) in the form of a 
40x40 matrix. Then, each sample is converted to the radiance as a string of 1600 
values. Images recorded with the original ifSOO system were reduced using somewhat 
smaller sample sizes to cover only sea surface. 

In the analysis of these quantitative sample values, the first set of 
calculations of basic statistics were performed by using the STATGRAPHICS 
progreun (Statgraphics is a registered trademark of Statistical Graphics 
Corporation). These calculations are 

* average 

* variance 

* standard deviation 

* geometric mean 

* minimum, maximum and range 

* standard error 

* standardized skewness and standardized kurtosis 

* median 

* mode 
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The steindardized coefficients show the deviations from the normeil distribution 
while the average, median, mode and geometric mean measure the central tendency 
and standard deviation, range and veiriamce measure spread of the data. [Ref. 11] 

After initial analysis, the samples were normalized and a Chi-square test 
was conducted The results are presented in the following sections while these 
results are discussed in the next Chapter. 

2. Analysis Of Sky Background Noise 
a. SW- Clear Sky 

A total 39 images of the SW clear sky have been analyzed. These 
images were taiken from the roof of Spanagel Hall at the Naval Postgraduate School, 
Monterey, Ca during August and September, 1989.The first 19 images were 
recorded with a randomly selected scarmer direction. The scanner direction during 
the recording of the last twenty images was selected in such a way to examine the 
effect of the scattered sunlight on the distribution type. This was done by choosing a 
portion of sk>’ with strong sunlight influence. 

A plot of the eimplitude distribution for all 39 images is shown in 
Figure 12. The three features seen in this figure eire due to differences in the mean 
amplitude This veiriation of the mean is primarily due to chemges in the azimuth 
and elevation of the viewing direction (see Table 4). As a result of the scatter in the 
data seen in this Figure, the mean value of each image was subtracted from the data 
producing normalized amplitudes (/i=0). The results aire shown in Figure 13. The 
full Normalized sample values (with /i = 0 eind (T =1) eire shown in figure 14. Figure 
15 shows the histogram of sample values from an image data and the standard 
Normal distribution. 
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The Chi-square test was conducted for each individual image as 

well as the overall data set under the null hypothesis that the underlying 

distribution was Normal. In order to test this the 12 cells with a width of 0.5 were 

chosen. The Chi-square values, observed and expected cell counts 2 ffe listed in 

Table 5 for the overall data set. The computed is 15.53218. The tabulated 

Chi-square for k — 12, eind m — *l, and a significance level of 0.05 is = 16.919. 

-05, fi 

(Recadl from Chapter 2, ib-1— m is the degrees of freedom, k is the number of cells, m 
is the number of pzu'ameters estimated.) Since 15.532 is less than 16.919, a normal 
distribution provides an excellent fit to the data. Observed euid expected cell counts 
are illustrated in Figure 16. Table 6 shows the anzilysis results of the SW clear sky. 




Figure 12. The Scatter Plot of Seunple Values of SW Cl. Sky (10"^ W cm"^ sr"') 
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Figure 13. The Scatter Plot of Szimple Values With Zero Mean (SW Clear Sky) 
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Figure 15. Standard Normal Eiistribution and Histogram of a Sample (SW CL. Sky) 
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Table 5. Cm-SQUARE TEST (SW CLEAR SKY) 



Cell 


Lower 


Upper 


Observed 


Expected 


Chi-square Test 


Num. 


Limit 


Limit 


Cell Count 


Cell Count 


Values 


1 


< 


-2.5 


9.2450 ■ 


9.9 


043336 


2 


-2.5 


-2.0 


22.9745 


26.5 


0.469008 


3 


-2.0 


-1.5 


61.8166 


70.5 


1.069507 


4 


-1.5 


-1.0 


132.7684 


147.0 


1.377809 


5 


-1.0 


-0.5 


253.8439 


239.8 


0.822478 


6 


-0.5 


0 


355.0097 


306.3 


7.7460999 


7 


0 


0.5 


290.3175 


306.3 


0.8339499 


8 


0.5 


1.0 


215.2746 


239.8 


2.508327 


9 


1.0 


1.5 


146.5009 


147.0 


0.001695 


10 


1.5 


2.0 


72.1412 


70.5 


0.038208 


11 


2.0 


2.5 


27.8728 


26.5 


0.071117 


12 


2.5 


< 


12.2348 


9.9 


0.550651 


h 


=12. /l: 


=2, k-i 


1-1=9 


Chi-square = 15.53218 




Figure 16. Observed and Expected Cell Counts (SW Clear Sky) 
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Table 6. THE RESULTS OF SW - CLEAR SKY ANALYSIS 



Average 


: 25.2716 x 10"J 


i W cm-2 sr-i 


Variance 


: 0.4 ( 


f/iW cm"^ sr"*)2 


Standard dev. 


: 0.63 


/iW cm-2 sr-* 


Minimum 


: 23.29 x 10‘S 


W cm"2 sr"* 


Maximum 


; 30.3 x 10-5 


W cm-2 sr-‘ 


Range 


: 7.01 x 10-5 


W cm"2 sr"‘ 


Average Range 


: 0.41 x 10-5 


W cm-2 sr-‘ 




Goodness of Fit 




Number of Samples : 39 




Sie. 


Level > 0.01 : 16 


41% 


Sig. Leel > 0.005 : 19 


48% 



6. SW - Cloudy Sky 

A total of 8 images of the SW cloudy sky were analyzed. The first 6 
images were recorded randomly from overcast sky with high altitude clouds. Two 
images were recorded from 80% cloudy sky with medium altitude clouds. The 
purpose of recording these last two images W2is to see the effect of the two different 
kinds of background noise sources on the amplitude distribution. The normzJized 
sample values from the 80% cloudy sky axe shown in Figure 17. A discussion of this 
figure is deferred to the next chapter. 

The normalized szunple values from the overcast sky are shown in 
Figure 18. The Chi— squared test was conducted under the same null hypothesis 
described in the previous section. The computed Chi-squaie is 17.34588. The 
tabulated Chi— squeire for i = 12 , m = 2 and a significance level of 0.025 is 
19.022. Since 17.34588 is less than 19.022, a normal distribution provides a good fit 
to data. The observed and expected cell counts are listed in Table 7 and illustrated 
in Figure 19. Table 8 shows the euialysis results of the SW samples of the overcast 
sky. 
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Figure 17. Normalized Sample Values of 80% Cloudy Sky (/i = 0 and <r = 1) 



100.00 



- 2.00 0.00 2.00 

Normalized Amplitude 



Figure 18. Normalized Sample Values of Overcast Sky (/i = 0 zind (7=1) 
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Table 7. CHI-SQUARE TEST (SW OVERCAST SKY) 



Cell 


Lower 


l^per 


Observed 


Expected 


Chi-square Test 


Num. 


Limit 


Limit 


Cell Count 


Cell Count 


Values 


1 


< 


-2.5 


4.b6t)7 


9.9 


27766542 


2 


-2.5 


- 2.0 


24.3333 


26.5 


0.177149 


3 


- 2.0 


-1.5 


82.6667 


70.5 


2.099685 


4 


-1.5 


- 1.0 


142.1667 


147.0 


0.158919 


5 


- 1.0 


-0.5 


238.5000 


239.8 


0.001046 


6 


-0.5 


0 


310.0000 


306.3 


0.052288 


7 


0 


0.5 


280.1667 


306.3 


2.180919 


8 


0.5 


1.0 


238.1667 


239^ 


0.002906 


9 


1.0 


1.5 


155.3333 


147.0 


0.472411 


10 


1.5 


2.0 


93.6667 


70.5 


7.612687 


11 


2.0 


2.5 


24.5000 


26.5 


0.150943 


12 


2.5 


< 


5.83333 


9.9 


1.670483 


ifc = 12 , /»; 


.2, k-i 


»-l=9 


Chi-square = 17.34588 
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Table 8. THE RESULTS OF SW - OVERCAST SKY ANALYSIS 



Avei^e 

Variance 

Standard dev. 

Minimum 

Maximum 

Range 

Average Range 



: 28.95 X lO’® 
: 1.75 

: 1.34 

: 25.40 X 10-s 
: 35.73 X lO'^ 
: 10.33 X 10-5 

: 0.78 X 10-5 

Goodness of Fit 



W cm"2 sr"* 
(/tW cm“2 sr"‘)2 
fiM cm" 2 sr"* 

W cm"2 sr"> 

W cm-2 sr-i 
W cm"2 sr"* 

W cm-2 sr-» 



Number of Samples : 6 

Sig. Level ^0.01 : 1 17X 

Sig. Level > 0.005 : 3 507. 



c. LW- Clear Sky 

A total of 30 images have been analyzed for the LW clear sky. 

Some of the images were recorded under a partly cloudy weather condition but the 

data were taken only from the clear part of the sky. The normalized sample values 

from these 30 images were shown in Figure 20. This normalized scatter plot 

indicates that the sample VEilues are piled about the meem. 

The result of the Chi— squre test showed that normal distribution 

did not provide a good fit to data under the null hypothesis described in the 

previous sections. The computed Chi-square is 174. The tabulated Chi— square for k 

= 12 , m = 2 and a significance level of 0.005 is <1^ = 23.587. Since 174 is much 

® - 005,9 

more than 23.587, the null hypothesis is rejected. The deviation of observed cell 
counts from expected cell counts can be seen easily in Figure 21. None of the 30 
samples individually fitted to normal distribution either. Table 9 shows the analysis 
results of the LW samples of the clear sky. The results of this section are discussed 
in the next chapter 



31 




Figure 20.The Plot of Normalised Sample Values (LW Cl. Sky /i = 0 and <r = 1) 
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Table 9. THE RESULTS OF LW - CLEAR SKY ANALYSIS 



Average 


: 22.55 X lO"* 


V cm"2 sr"^ 


Variance 


: 1.21 X 10-10 


(W oq-2 sr'O^ 


Standard dev. 


: 0.11 X 10--* 


V cm"2 sr"^ 


MiniflUffl 


: 22.55 X lO-* 


V cm-2 sr-» 


Maxinniis 


: 35.05 X lO--* 


V cm-2 sr"^ 


Range 


; 12.50 X 10-< 


V ca-2 sr"^ 


Average Range 


: 0.49 X lO--* 


V cm-2 sr"^ 




Goodness of Fit 




Number of Samples : 30 




Sig. Level > 0.01 : NONE 


Sig. 


Level > 0.005 : NCB® 



3. Aneilysis Of Sea Background Noise 
a. SW — Sea Surface 

A total of 20 images have been analyzed for the SW sea surface. 
The images were teiken from the Hopkins Marine Station on Sept. 14, 1989 Due to 
hmitations in the field— of— view’, the images included objects other than the sea 
surface In order to cover only the sea surface, eight hundred pixels of image data 
w’ere selected The normalized sample values from these 20 images w’ere shown in 
Figure 22. 

The result of the Chi-square test show’ed that normal distribution 

did not provide a good fit to data for overall sample values. The computed 

Chi— square is 70.88 The tabulated Chi— square for ^ = 12 , m = 2 zmd a 

sigmficeince level of 0.005 is = 23.587 Since 70.88 is much more than 23.587, 
® -OOSiS 

the null hypothesis is rejected. The de\uation of observed cell counts from expected 
cell counts can be seen in Figure 23. However, four samples fitted to a normal 
distribution individually Table 10 shows the analysis results of the SW samples of 
the sea surface 
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Figuie 23. Observed and Expected Cell Coxmts (SW Sea Surface) 
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Table 10. THE RESULTS OF SW - SEA SURFACE ANALYSIS 



Average 

Variance 

Standard dev. 

Mininum 

Maximum 

Range 

Average Range 



: 26.45 X lO"* 

; 5.27 X lO-is 

: 7.26 X 10-7 

: 25.40 X 10-5 
: 28.40 X 10-5 
: 3.01 X 10-5 

: 0.55 X 10-5 

Goodness of Fit 



V sr-^ 

(W cm -2 8 r -»)2 
W cm-2 sr-‘ 

W cffl-2 sr-i 

V cffl-7 sr-^ 

W cm-2 sr-* 

V sr-^ 



Number of Saiq)le8 : 20 

Sig. Level > 0.01 : 4 20Z 

Sig. Level > 0.005 : 5 257. 



b. LW — Sea Surface 

A toted of 20 images of the LW sea surface were analyzed. 

The images were taken during November, 1988, from the Moss Landing Meirine Lab 

[Ref 2] The sample size was reduced to 1279 in order to cover only sea surface. 

The normalized sample values from these 20 images were shown in Figure 24. 

Chi-square values, eind observed emd estimated cell counts are 

listed in Table 11. The computed Chi-squzire is 22.75541. The tabulated 

Chi-square for A = 12, m = 2, and a significemce level of 0.025 is = 23.584. 

■ 005>9 

Since 22 75 is less than 23.584, a normal distribution provides a good fit to the data. 
Observed and expected cell counts Eire also illustrated in Figure 25. Table 12 shows 
the analysis results of the LW samples. 
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Table 11. CHI-SQUARE TEST (LW SEA SURFACE) 



Cell 


Lower 


I^per 


Observed 


Expected 


Chi-square Test 


Num. 


Limit 


Limit 


Cell Count 


Cell Count 


Values 


1 


< 


-2.5 


4.76 


7.9 


06D3S? 


2 


-2.5 


-2.0 


11.95 


21.17 


4.013383 


3 


-2.0 


-1.5 


62.65 


56.3 


0.713381 


4 


-1.5 


-1.0 


132.45 


117.4 


1.928489 


5 


-1.0 


-0.5 


211.60 


190.9 


2.244293 


6 


-0.5 


0 


215.20 


244.65 


3.546687 


7 


0 


0.5 


267.05 


244.65 


2.049566 


8 


0.5 


1.0 


171.05 


190.9 


2.191969 


9 


1.0 


1.5 


103.60 


117.4 


1.625744 


10 


1.5 


2.0 


62.45 


56.3 


0.66907 


11 


2.0 


2.5 


24.25 


21.17 


0.449082 


12 


2.5 


< 


11.95 


7.9 


2.066646 


fc=12, n: 


=2, k-i 


»-l=9 


Chi-square = 22.75541 




Figure 24.The Plot of Normalised Sample Values (LW Sea Surf., /i = 0 and <r = 1) 
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Table 12. THE RESULTS OF LW - SEA SURFACE ANALYSIS 



Average 


: 43.00 X lO--* 


V cm-* sr-i 


Variance 


: 3.14 X 10-1* 


(V cm-* sr-*)* 


Standard dev. 


: 0.56 X 10-5 


w cm-* sr-i 


Minimnn 


: 42.75 X 10-* 


V cm-* sr-i 


Maximum 


: 43.30 X lO-i 


V cm-* sr-i 


Range 


: 0.55 X lO--* 


V cm-* sr-i 


Average Range 


: 0.33 X 10 -* 


V cm-* sr-i 




Goodness of Fit 




Number of Samples : 20 




Sig. Level > 0.01 ; 6 


307, 


Sig. Level > 0.005 : 10 


50X 




Figure 25. Observed and Expected Cell Counts (LW Sea Surface) 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The results presented in chapter IV show that the amplitude distributions of 
several infrared backgrounds in the spectral region 3—5.6 /im and 8—14 fan obey a 
normal distribution. In the 3—5.6 fan region, a normal distribution provides a good 
fit at the 0.05 level of significance for the clear sky background. The percentage of 
individual images obe)hng a normal distribution at the 0.01 level of significance is 
41% and 48% at the 0.005 level of significemce. The data show the mean amplitude 
varies with the viewing direction (Figure 12). However, the variemce of the SW clear 
sky images does not change much. This cein be seen by comparing Figures 13 eind 
14. 

In the same spectral region, the overall data set for the overcast sky also fit a 
normal distribution at the 0.025 level of significance. One of the six images fitted a 
normal distribution at the 0.01 level of significance while three images fitted a 
normal distribution at the 0.005 level of significance. However, the distribution of 
the images containing both cloud and blue sky appears to deviate from a Gaussian 
[Fig. 17]. It is believed that this deviation is the result of a superposition of the two 
parent distributions. Based on the average amplitudes of the clear sky zind overcast 
sky (see tables 6 zind 8), it appears that the large peak in Fig. 17 is due to the 
overcast sky zind "shoulder" on the left is due the clear sky. 

In the 3—5.6 fan spectral region a normeil distribution does not provide a good 
fit to the overall data set for the sea background. Due to limitations in the 
field— of— view, the images included objects other than the sea surfaw:e. This may 
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cause a deviation in the sample values from the normal distribution. Nevertheless, 
20% of the sea surface images fit a normal distribution at the 0.01 level of 
significance and 25% of images fit a normal distribution at the 0.005 level of 
sigmficeince. This is cleeirly reflected in the overall amplitude distribution shown in 
Figure 23. 

The aunplitude distribution of the sea backgroimd in the spectrad region 
8—14 /zm obeys a Gaussion distribution at the 0.005 level of significance. The 
percentage of individuad images obeying a normad distribution is about 30% at the 
0.01 level of significance and is about 50% at the 0.005 level of significance. The 
clear sky background in the spectral region 8—14 fan does not appear to obey a 
Gaussian distribution. The data collected from the clear sky, however, may not be 
valid because the amplitudes of the saimple values from this band are very close to 
the lower end of the dynamic range of AGA THV 780 scanner. Due to this possible 
non-linearity, the measured sample values of LW— clear sky may be distorted. 

The results of this investigation are summarized in Table 12. The amplitude 
distribution of the clear sky (3—5.6 ^m), overcast sky (3—5,6 fan) and sea surface 
(8—14 fan) backgrounds can be described by Gaussian distribution. No conclusion is 
made regarding the parent distribution for the clear sky (8—14 fan) or the sea 
surfeice (3—5.6 fan) beickgrounds. Where a sample set extends over two background 
types (eg. cloud and clear sky) the distribution appears to be a bimodaJ one 
consisting of a superposition of the separate distributions for the two types. 
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Table 12. SUMMARY OF THE RESULTS 



Background 


3-5.6 /xm 
Variance* 


Distribution^ 


8-14 /xm 
Variance 


Distribution 


Clear 

Sky 

Overcast 

Sky 

Sea 

Surface 


0.4 X 10’*2 
1.75 X 10-*2 
0.53 X 10-*2 


Normal (0.05) 

Normal (0.025) 

Withhold 

Judgement 


1.21 X 10’ *0 


withhold 

Judgement 


3.14 X 10-‘» 


Normal(0.005) 



1 Unit; (W cm-2 81-1)2. 

2 Value in parenthesis indicates significance level. 



B. RECOMMENDATIONS 

Although the parent distributions of several important infrared beu:kground8 
have been modeled with a Gaussian distribution, estimation of the mean and the 
standard deviation must be done very carefully to get the output of this modeling 
close to the real life vedues. In order to do this, the relationship between these two 
parameters and various background pareimeters, such as the air and the sea 
temperature, should be known. In general, the field— of— view will be composed of 
two or more background types. In this case, it is necessary to know the type and 
magnitude of the individual noise sources within a composite background. This 
requires further statistical analysis of the IR background noise parameters assuming 
the amplitude distribution is Gaussian. 
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APPENDIX A 



C 

C 

C 

C 

C 

C 

C 
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LISTING OF THE SBl FORTRAN PROGRAM 
**-***-^-*PURPOSE’****’*’^* 

THIS PROGRAM CALCULATES IN-BAND FLUX. 

INPUTS; LOWER AND UPPER LIMITS OF BAND 

THE THERMAL VALUE (lU), OR TEMPERATURE IN K 
OR OC 

NUMBER OF INTERVALS 
OUTPUT W/cm^, Photon/cm^ vs lU 
INTEGER N.M, I, J,K,S, RET, BAL, DAL, YES, L.BIL 
REAL LAM,T(1000),ARA,A,G,C,Q(1000),U(1000),C1,C2,IX, 
+V4.W(1000),C3,LAM1,LAM2,DIV,TMAX,TMIN,TRAN,V1,V2,V3 
OPEN(UNIT=l, FILE='SONUC.TXT', 

+STATUS='UNKNOWN', ACCESS='SEQUENTIAL') 

Cl = 37415 
C2 = 14387.9 



C3 = 188365 
IX = 1000 



IX = IX^-6 
C3 = IX*C3 

PRINT", ' ENTER WAVELENGTH 1 (lOxMICRONS) ' 
READ ", LAMl 

PRINT *, ’ ENTER WAVELENGTH2 (lOxMICRONS) ' 
READ ", LAM2 



41 



PRINT *, 'ENTER NUMBER OF INTERVALS OF THE INTEGRATION' 
READ DIV 

ARA = (LAM2-LAM1)/DIV 
WRITE(*,222) ara 
222 FORMAT(3X,G15.8) 

IF(LAM1.LT.7) THEN 

A=811203 

G=31 17.74 

C=1 

GOTO 2 
ENDIF 
A = -3581 
G = 1506.49 
C = -0.436 

2 PRINT *, 'DO YOU USE TMV 780 UNITS' 

PRINT *, ' ENTER 1 FOR YES ' 

PRINT ' ENTER 0 FOR NO ' 

READ ♦. YES 

IF(YES.GT.O) GOTO 3 

PRINT *, ' WHAT IS THE INPUT UNIT? ' 

PRINT *, ' ENTER 1 FOR TEMPERATURE IN KELVIN ' 

PRINT ' ENTER 0 FOR TEMPERATURE IN CELCIUS' 

READ *, DAL 
IF(DAL.LT.l) THEN 

PRINT *, ' ENTER MIN TEMPERATURE IN CELCIUS' 

READ TMIN 
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TMIN = TMIN+273 

PRINT *, ' ENTER MAX TEMPERATURE IN CELSIUS' 

READ *, TMAX 

TMAX = TMAX+273 

GOTO 1 

ENDIF 

PRINT *. ' ENTER MIN TEMPERATURE IN KELVIN' 
READ *, TMIN 
TMIN = TMIN 

PRINT ' ENTER MAX TEMPERATURE IN KELVIN' 
READ ^ TMAX 
TMAX = TMAX 
GOTO 1 

C 

3 PRINT =‘, ' ENTER MIN ISOTHERM UNIT ' 

READ *, TMIN 
TMIN = TMIN 

PRINT ^ ' ENTER MAX ISOTHERM UNIT ' 

READ *. TMAX 
TMAX = TMAX 

C 

1 TRAN = TMAX-TMIN 
K = INT(TRAN) 

K=K+1 
DIV = DIV+1 

C 
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DO 600 J=1,K 
W(J) = 0 

IF(YES.GT.O) THEN 

U(J) = TMIN+J-1 

T(J)=(G/(LOG(((A/U(J))+l)/C))) 

GOTO 4 
ENDIF 

T(J) = TMIN+J-1 
4 LAM = LAMl 

Q^itcJtc**************************************:*^********* 

DO 601 1=1, DIV 

VI = (C1/(LAM*LAM*LAM*LAM*LAM)) 

V2 = (C2/(LAM*T(J))) 

V3 = 1/(EXP(V2)-1) 

V4 = (V1*V3)*ARA 
LAM = LAM+ARA 
IF(I.EQ DIV) THEN 
V4 = (V4/2) 

ENDIF 

W(J) = V4+W(J) 

601 CONTINUE 
600 CONTINUE 

no D0 700J=1,K 

Q(J) = 0 

IF(YES.GT.O) THEN 
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U(J) = TMIN+J-1 
T(J)=(G/(L0G(((A/U(J))+1)/C))) 

GO TO 5 
ENDIF 

T(J) = TMIN+J-1 
5 LAM = LAMl 

DO 701 1=1, DIV 

VI = (C3/(LAM’^LAM*LAM*LAM)) 

V2 = (C2/(LAM*T(J))) 

V3= 1/(EXP(V2)-1) 

V4 = (V1*V3)*ARA 
LAM = LAM+ARA 
IF(I EQ.DIV) THEN 
V4 = (V4/2) 

ENDIF 

Q(J) = V4+Q(J) 

701 CONTINUE 
700 CONTINUE 
111 J=1 

K=K+1 

300 IF(J.LT K) THEN 

WRJTE(1,103) U(J),T(J),Q(J),W(J) 

WRJTE(M03) U(J),T(J),Q(J),W(J) 

J=J+1 

103 F0RMAT(3X,1F4 1,1X,1F10.3,1X,1G12.4,1X,1F15.8) 
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GO TO 300 
ENDIF 

PRINT ♦, 'DO YOU WANT TO RUN THE PROGRAM AGAIN' 
PRINT 'ENTER 1 FOR YES ' 

PRINT 'ENTER 0 FOR NO ' 

READ *, BIL 
IF(BIL.LT.l) GOTO 998 
GO TO 997 
998 STOP 
END 
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APPENDIX B 



LISTING OF THE SB2 FORTRAN PROGRAM 
C ******’*PURPOSE****’*** 

C THIS PROGRAM READS THE TIC-8000 IMAGE FILE AND 
C COMPUTES THE VALUE OF EACH PIXEL AS ^C, W/cm2 or lU. 

C INPUT IS READ FROM THE FILE CALLED INPUT.IMG 
C SPOT VALUES CAN BE READ FROM SCREEN 
C OUTPUT CAN BE WRITTEN IN A FILE CALLED SB20UT.MAT 
C AS A 40x40 MATRIX 

INTEGER N,M, I, J,K,S, RET, BAL, DAL, YES, L, IX, BIL 
INTEGERS B(20446) 

REAL RAN,D(140,140),LEV,A,G,C,W,0 
998 0PEN(UNIT=99,FILE='a\input.IMG', 

+ACCESS='DIRECT',STATUS='UNKNOWN',FORM='BINARY', 

+RECL=20446) 

OPEN(UNIT=l , FILE='SB20UT.MAT’ , STATUS='UNKNOWN' ) 
PRINT ^ ' HOW DO YOU WANT OUTPUT ' 

PRINT *, ' ENTER 0 FOR ISOTHERM UNITS ' 

PRINT *, ' ENTER 1 FOR TEMPERATURE (OC) ' 

PRINT ^ ' ENTER 2 FOR IRRADIANCE (W/cm2) ' 

READ *, BAL 
IF(BAL EQ.O) THEN 
GO TO no 
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ENDIF 

PRINT ♦, ' ENTER 0 FOR SW ' 

PRINT ♦, ’ ENTER 1 FOR LW 
READ DAL 
IF(DAL.LT.1)THEN 
A=811203 
0=3117.74 
0=1 
W=4.69 
W=W*(.00001) 

0=1 

GOTO no 
ENDIF 
A = -3581 
G = 1506.49 
C = -0.436 
W = 3.03 
W=W/(10000) 

O = 3.624 

no WRITE(*,222) W 
222 FORMAT(3X,G15.8) 

PRINT*,' ' 

PRINT*,' ' 

PRINT *, ' READING IMAGE FILE ' 
READ (99,REC=1) (B(I), 1=1,20446) 

102 FORMAT(19600A1) 
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PRINT ' ENTER RANGE ' 

READ *. RAN 

PRINT ' ENTER LEVEL ' 

READ *, LEV 
PRINT ' 

PRINT 'PLEASE WAIT ' 

1=846 

DO 601 J=l,140 
DO 602 K=l,140 
1 = 1+1 

IF(B(I).LT.O) THEN 
D(J,K)=FLOAT(B(I)) 

D(J,Kj=256+D(J,K) 

D(J,K)= (((RAN"((D(J,K)-128.5)))/255)+LEV) 

IF(D(J,K).LT.O) THEN 

D(J,K) = 0.000001 

ENDIF 

GO TO 111 

ENDIF 

D(J,K)=FLOAT(B(I)) 

D(J,K)= (((RAN’‘((D(J,K)-128,5)))/255)+LEV) 
IF(D(J,K).LT.O) THEN 
D(J,K) = 0.000001 
ENDIF 

111 IF(BAL.EQ 1) THEN 

D(J,K)=(G/(LOG(((A/D(J,K))+l)/C))) 
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D(J,K)=(D(J,K)-273) 

ENDIF 

IF(BAL.EQ,2) THEN 
D(J,K)= W*(D(J,K)+0) 

ENDIF 

602 CONTINUE 
601 CONTINUE 

113 PRINT *, ' DO YOU WANT SPOT VALUE ' 
PRINT *, ' ENTER; YES = 1, NO = 0 ' 

READ YES 
IF(YES.GT.O) THEN 

PRINT *, ' THE COORDINATES OF THE PIXEL' 
PRINT ' ' 

PRINT ' ENTER X, 0<X<139 ' 

READ*, IX 
PRINT *. ' ' 

PRINT *, ' ENTER Y, 0<Y<139, (From Down) ' 

READ *, M 

S=M 

M=(140-M) 

L=M/2 

L=2*L 

IF(M.GT.L) THEN 
J=((M/2)+71) 

GO TO 114 
ENDIF 
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114 



J=(M/2) 

K=IX+1 



WRITE(*,115) IX,S,D(J,K) 

115 F0RMAT(1X, '('.113.3, ',',113.3,') =',F15.8) 

GOTO 113 
ENDIF 

PRINT *, ' DO YOU WANT OUTPUT FILE ' 

PRINT *, ' ENTER; YES = 1, NO = 0 ' 

READ NO 
IF(x\O.LT l) GOTO 999 

117 IF(BAL.GT.O) THEN 

109 PRINT *, 'OUTPUT WILL BE WRITTEN IN SB20UT.MAT' 
PRINT * 

PRINT " 

PRINT * 

PRINT ■ 

PRINT * 

PRINT ’ 

DO 116 5=1,100000 

118 CONTINUE 
GOTO 116 
ENDIF 

PRINT 'OUTPUT WILL BE WRITTEN IN SB20UT.MAT' 
PRINT r ' ' 



PRINT ' 
PRINT - ' ' 
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PRINT *, ' ' 

PRINT *, ' ' 

DO 119 8=1,100000 
119 CONTINUE 
116 M=26 

300 IF(M.LT.46) THEN 

J=M+70 

WRITE(1,103) (D(J,K),K=51,90) 

WRITE(*,103) (D(J,K),K=51,90) 

J=M 

WRITE(1,103) (D(J,K),K=51,90) 

WRJTE(*,103) (D(J,K),K=51,90) 

M=M+1 

103 FORMAT(140F15.8) 

GO TO 300 
ENDIF 

999 CLOSE(l) 

PRINT ’ DO YOU WANT TO RERUN THE PROGRAM ' 
PRINT *, ’ ENTER; YES = 1. NO = 0 ’ 

READ YES 
IF(YES.LT.l) GOTO 997 
CLOSE(99) 

CLOSE(l) 

GOTO 998 
997 STOP 
END 
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APPENDIX C 



IMAGE TRANSFER PROCEDURE 

A. GENERAL 

The thermeJ image data recorded and stored on the diskettes using the ifSOO 
computer were transferred to the IBM— AT Computer using the DISPORT program 
via a RS-232 interface cable (DISPORT is a registered trademark of GESOTEC 
Soft— cind Heirdwzire GmbH, FRG). The transmission was ceirried out through the 
serieil ports of the computers using ASCII format. PROCOMM PLUS softweire 
program was used to send required commeinds to the ifSOO computer (PROCOMM 
PLUS is a registered trademark of DATASTORM TECHNOLOGIES, INC.). 
However, any communication program can be used for this purpose. 

Data transfer is carried out in blocks. There are two types of blocks. The first 
type is a block of 252 ASCII characters for the image directory. The organization of 
the image directory is shown in Table 13. The second type is the thermal image 
information. The 128 columns by 64 rows of the thermal image are transferred 
sequentially as bytes starting at the upper left-hand corner of the image The 
thermeil information of a image is broken down into 64 blocks. Eaw:h block has 384 
ASCII characters. A code containing three ASCII characters represents the value of 
one byte. The most significant digit is transferred first and the least significzmt digit 
is sent last. Thus, each block contains 128 bytes which represent one row of the 
thermal image. [Ref. 12] 

For a complete description see the Disport Operating Manual [Ref. 12]. 
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TABLE 13. ORGANIZATION OF THE IMAGE DIRECTORY 
(Reproduced from Ref. 7) 















Param. I 


Char. No . 


I stored 


as : 


I comment 


disk I 

label I 


1.2 


I Ascii 
I 




I 

I 


im.no. I 


3-8 


I 16 


bit 


int 


I image number 


RANGE I 


9-14 


I 16 


bit 


int 


I 


LEVEL I 


15-20 


I 16 


b1 1 


int 


I 


abs(A) I 


21-29 


I 24 


bit 


int 


I abs(A)=65536*DI(5)+256*DI(6)+DI(7) 


B I 


30-35 


I 16 


bit 


int 


I 


abs(C* I 
1000) I 


36-44 


I 24 
I 


bi t 


int 


I abs(C)= (65536*01 ( 9 )+256*DI (10)+ 
I +DI (11) )/1000 


Dcal*10 I 
I 


45-50 


I 16 
I 


bi t 


int 


I calibration distance (m) 
I dcal=0I (9)/10 


sign I 

(A.C)I 


51-56 


I 16 
I 


bit 


int 


I bit 0: sign of C 
I bit 1 : sign of A 

j 


EPS*100 I 


57-62 


I 16 


bit 


int 


I emissivity obj. * 100 


DIS 1 


63-68 


I 16 


bit 


int 


I object distance (m) 


Tamb*10 I 


69-74 


I 16 


bit 


int 


I ambient temperature (K) * 10 


Tatm*10 I 


75-80 


I 16 


bi t 


int 


I atmospheric temperature (K) * 10 


TRMS*100I 


81-86 


I 16 


bit 


int 


I transmission * 100 


alpha* I 
1000001 


87-95 


I 24 
I 


bit 


int 


I alpha=( 65536*01 (19)+256*DI (20)+ 
I +01(21))/100000 


beta* I 

1000001 


96-104 


I 24. bit 
I 


int 


I beta=(65536*OI(22)+256*OI(23)+ 
I +0I(24))/100000 

I 


0/R sw I 
I 


105-110 


I 16 

I '■ 


bit 


i nt 


I direct/relative switch 
I 1 : rel ati ve, 0 : di rect 


EPSR*100I 


111-116 


I 16 


bi t 


i nt 


I emissivity ref., epsr=0I ( 14 ) /lOO 


XR I 


117-122 


I 16 


bit 


int 


I horizontal position ref. 


YR I 


123-128 


I 16 


bit 


int , 


I vertical position ref. 


VALR I 


129-134 


I 16 


bit 


int 


I digital value of reference point 


Tref*10 I 

j 


135-140 


I 16 

T 


bit 


int 


I reference temperature (K) * 10 

j 


SCTY I 


141-146 


I 16 


bit 


int 


I scanner type 


SEVER I 


147-149 


I 


ascii 


I scanner version 


SERNO I 


150-155 


I 16 


bit 


int 


I serial no. 


FILCOD I 


156-158 


I 


asci i 


I filter code 


F0V*2 I 


159-164 


I 16 


bit 


int 


I fov=0I(39)/2 


APERT*10I 


165-170 


I 16 


bi t 


int 


I aperture, ap=0I(40)/10 


RN0*10 I 


171-176 


I 16 


bit 


int 


I revision no., rno=0I ( 41 ) /lO 


FI/FRSW I 
I 


177-182 


I 16 
I 


bi t 


int 


I FIELO/FRAME switch (l=Field, 
I 2=Frame) 


image I 

label I 


183-252 


I ascii 
I 




I 

I 
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B. TRANSFER PROCEDURE 

The following is a step-by-step description of the transfer procedure. 
Figure 26 shows the transfer order of a image data. 



* Connect the serial ports of the computers with a RS232 interface cable. 
There is no need to switch hnes in this cable. A straight pin— to-pin 
connection is all that necessary. 

Turn on both computers. 

* Configure the serial port of the ifSOO computer by \ising the 
CONFIG. The standard configuration is 

Beuid Rate: 9600, Data Bits: 8, Stop Bits: 1, Parity: NO. 

* Load the communication progreim to the IBM— AT computer and configure 
the serieJ port as described in the previous step. 

* According to the user manuaJ of the communication program, open a file 
to write the received ASCII characters. 

* Load DISPORT program to the ifSOO computer. After a short user dialogue, 
if800 is ready to send. 

* In the IBM— AT, be sure you are on the communication screen and everything 
written onto the screen will be sent to the ifSOO and will be the content of the 
file opened. From now on, do not to make any mistake beceoise each key stroke 
win be written into the file as zm ASCII character and will mix with the 
received image data ASCII characters. 

* Send an STX command to the ifSOO computer from IBM— AT computer. This 
can be done by typing either Ctrl B or Alt 002 (from numeric keyboard). 

* ifSOO responds and begins sending the image directory (252 ASCII characters). 
See Table 13. At the end of this block the ifSOO sends an ETX command 
(ASCII value: 003). 

* Send the second STX command. IfSOO begins to send the thermeil image 
information (384 ASCII cheiracters per block). At the end of each block the 
ifSOO sends an ETX command and waits for the next STX command. 

* Send an STX commeuid adter each ETX command. This procedure continues 
until 64 blocks of the thermal image have been transferred. 

* At the end of the imaige trainsfer, DISPORT terminates and ask for another 
data transfer. 

* In the IBM— AT close and save the file. If you wauit to trainsfer ainother image 
data, open a new file emd repeat the transfer steps. 
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